Abstract. On the basis of multislit spectroscopy of 180 
stars in the ionising cluster of 30 Doradus we present re- 
liable radial velocities for 55 stars. We calculate a radial 
velocity dispersion of ~ 35 km s _1 for the cluster and 
we analyse the possible influence of spectroscopic binaries 
in this rather large velocity dispersion. We use numerical 
simulations to show that the observations are consistent 
with the hypothesis that all the stars in the cluster are bi- 
naries, and the total mass of the cluster is ~ 5 x 10 5 M . 
A simple test shows only marginal evidence for dynamical 

. mass segregation which if present is most likely not due 

' to dynamical relaxation. 

Key words: Stars: early-type - stars: kinematics - bina- 
O , ries: spectroscopic - galaxies: clusters: general - galaxies: 
' Magellanic Clouds. 
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1. Introduction 

A significant fraction of the old stars we now observe in 
galaxies appear to have formed in Starbursts. Therefore, 
understanding violent star formation becomes crucial if 
we want to understand the star formation history of the 
Universe. 30 Doradus in the LMC is the nearest and best 
studied example of a massive starburst cluster and thus it 
has become a sort of 'Rosetta Stone' for deciphering the 



physics of violent star formation processes (Walborn 1991 
Selman et al. 1999b| (Paper III) and references therein). 



Although 30 Dor has been the subject of intensive ob- 
servational effort from the ground and space there are still 
a number of critical problems that remain unsolved. Per- 
haps the most burning open problem is the iss ue of mass 
segregation first raised by Malumuth & Heap ( 1995 ) that 
has important implications for our understanding of the 
process of cluster formation in general (Clarke 2001). 

This Paper is the fourth of a series devoted to a 
comprehensive study of the 30 Doradus starburst cluster 
(NGC 2070). We can summarize the central results of the 
previous papers of this series as follows: 

— The Initial Mass Function (IMF) of NGC 2070 between 
3 and 120 M Q is very well represented by a Salpeter 
power-law. The total 'photometric' mass of the cluster 
- inside 20 pc - from 0.1 M to 120 M is M phot ~ 
3 x 10 5 M Q . 

— The star-formation history of the cluster shows three 
distinct peaks centered 5 Myrs, 2.5 Myrs, and < 
1.5Myrs ago with no evidence for a spatial segregation 
of ages. 

— The density distribution of stars between 20 and 
60 M Q is a power-law of slope of —2.85. This means 
that the total mass increases very slowly with radius 
(r 015 or even logarithmically. 
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— There is some evidence for mass segregation revealed 
by the fact that the most massive stars appear to be 
more concentrated toward the center. However, there 
is a distinct 'ring' of massive stars located at about 
6pc from the center of the cluster (R136). This ring 
contains stars of the three age groups described above. 

The purpose of the present Paper is to study the stel- 
lar kinematics of the brightest stars in the cluster using 
the NTT obse rvatio ns described in Paper II of this series 
(Bosch et al. 1999 ). We use the data to investigate the 
presence of dynamical mass segregation in the cluster, and 
to attempt to constrain the lower end of the IMF from the 
comparison between the photometric and the dynamical 
masses. 



2. Observations and Data Reduction 

The observations - already described in Paper II - were 
obtained with the La Silla NTT telescope using the Multi- 
slit option of EMMI's RILD mode. In this mode, grism # 
5 was used, which in combination with a 0.8 arcsec wide 
slit yields a resolution of 1.3 A / pixel. The wavelength 
range is 3600 to 6000 A in most cases, although the limits 
vary slightly between individual spectra, as they depend 
on the position of each slitlet within the mask. A total of 
7 masks were produced each including an average of 35 
slitlets. The basic reduction steps were described in Pa- 
per II including bias correction, flat fielding, wavelength 
calibration, and aperture extraction. Here we need to dis- 
cuss some aspects of the wavelength calibration not in- 
cluded in Paper II that are particularly relevant for the 
determination of radial velocities. 

During the observations a technical problem made the 
Argon lamp too weak to produce reliable identification 
lines for 2 masks while the strong He I 5875 A arc line 
is saturated for several spectra. The problem is worst when 
both effects appear at the same time as there are no reli- 
able comparison lines redwards Hei 5015 A. For example, 
this makes useless the Hen 5411 A absorption line for 
radial velocity measurements in the affected spectra. 
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The problems with the calibration lamps also altered 
the normal sequence in which the lamps were taken; the 
calibration exposures for the last three masks had to 
be taken at the end of the night. Although there were no 
telescope presets, the telescope was still tracking so the 
position angle of the instrument rotator slowly changed 
and this may introduce systematic errors. We address this 
issue in the following section. 

3. Radial Velocities 

3.1. Zero-point errors 

The problems with the calibration lamps may introduce 
systematic differences between the zero points of each 
mask. We have used the radial velocity of the [O i] 5577 
A Auroral line, whenever possible, to test for this effect. 

The straightforward approach would be to add alge- 
braically the radial velocity of the sky line to the veloc- 
ity of each star. Unfortunately, however, the prob- 
lems with the calibration lamps described above 
imply that several spectra have unreliable wave- 
length calibration in the region of [Oi]. 

For most spectra the formal (measurement) error in 
the position of the sky line is ~6,6 km s _1 . The mean 
radial velocity of the [O i] lines with reliable calibrations is 
-6.5 km s" 1 with a dispersion of 12.9 km s _1 . Thus, while 
we cannot correct the individual velocities for systematic 
zero point errors, the radial velocity dispersion must be 
corrected by this effect by subtracting quadratically the 
[Oi] dispersion a [Ql] = a/12.9 2 - 6.6 2 = 10.9 km s _1 . 

3.2. Velocity determination 

As mentioned above, the resolution of our combination 
of dispersion grating and camera yields 1.3 A per pixel. 
For 2-pixel sampling this corresponds to a spectral reso- 
lution of about 165 km s -1 at 4750A. The ultimate limit 
attainable in the precision of Doppler shifts is dom inated 
by the photon noise in the spectrum (Brown 199C ). The 
uncertainty in the measured radial velocity for the case of 
a single line of width w and depth d, measured in units of 

the continuum intensity 7 C , is given by 
cw 

dVms ~ \d(N w i c y/i (ij 

where N w is the number of samples across the width of the 
line, c is the speed of light, and A is the central wavelength 
of the line. For typical lines measured in our stars (A = 
4500A, u>=7A, N w =6, d=0.2, and I c = 2x 10 4 ), we obtain 
fo rms = 15 km s _1 . This is smaller than the velocity 
dispersion expected if the cluster is virialised and the total 
mass is close to our photometric estimate. However, the 
signatures of expanding stellar atmospheres and binaries 
may be much stronger than the virial motions. Thus, it is 
very important to constrain these effects with the data at 
hand before embarking in a high spectral resolution survey 
of the kinematics of the cluster. 



At our resolution the Balmer lines cannot be used to 
measure radial velocities because they are severely con- 
taminated by nebular emission. Therefore, we have re- 
stricted our analysis to stars with well exposed Hel and 
Hell absorption lines. In order to have an indication of the 
presence of atmospheric motions we have only considered 
stars with at least three He lines detected. This further 
restricts our sample to 97 spectra, several of which corre- 
spond to the same star. 

The centroids of the lines were determined from Gaus- 
sian fits using the package ngauss within IRAF. The fit- 
ting errors were used as weights to calculate the (weighted) 
mean velocity of each star. A conservative k — a filter was 
used to reject stars with suspected internal (atmospheric) 
motions. Thus, all stars with a dispersion of more than 
25 km s -1 between the measured lines were rejected. The 
final list is presented in Table [j] that gives, for each star, 
the Parker number, the spectral type from Paper II, posi- 
tion in arc-sec from the cluster center, assumed to be R136 



(Sclman et al. 1999a ), the weighted average radial velocity, 
and the weighted error. A number of stars appear to 
be binaries on the basis of showing double peaked 
lines (stars # 1024, 1369 and 1938), asymmetric 
line profiles (# 222, 613, and 1191), or different ra- 
dial velocities for the Hel and Hen lines (# 1998). 
These stars tend to have larger internal errors as can be 
seen in the second part of Table |l|. 
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Fig. 1. Spatial distribution of stellar radial velocities. Two 
stars are missing, as they fall out of the photometry area 
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Table 1. Stellar radial velocities 



Star id. 


Sp.Type 


M") 


AS(") 


(Vr) 


(Tint 


15 


08.5 V 


-107 


107 


234.9 


14.0 


32 


09 IV 


-102 


72 


272.1 


09.5 


124 


08.5 V 


-76 


25 


287.5 


06.9 


316 


06.5 V 


-50 


-164 


265.9 


09.5 


541 


07.5 V 


-29 


-65 


252.6 


07.1 


649 


08-9 V 


-20 


-106 


323.7 


09.4 


713 


05 V 


-15 


-53 


308.7 


11.8 


747 


06-8 V 


-13 


-142 


364.3 


22.8 


791 


05 V 


-09 


141 


310.7 


08.8 


805 


05-6 V 


-08 


-38 


292.1 


09.4 


863 


06.5 V 


-04 


-03 


308.0 


06.5 


871 


04 V((f )) 


-04 


-44 


290.3 


06.5 


901 


03 V((f )) 


-02 


26 


276.2 


08.7 


905 


O9-B0 V 


-02 


61 


198.2 


16.3 


975 


06-7 V((f)) 


02 


-27 


325.5 


05.9 


1022 


05: V 


04 


-139 


320.6 


07.2 


1063 


06-7 V 


06 


108 


267.0 


16.5 


1109 


09 V 


09 


-167 


229.7 


05.9 


1139 


B0 V 


11 


36 


225.5 


09.5 


1163 


04 If: 


12 


-72 


274.1 


11.3 


1247 


B0. 5 IV 


17 


91 


333.0 


11.6 


1283 


06 V:((f )) 


19 


-09 


352.2 


06.7 


1339 


B0-0.2 IV 


23 


-60 


265.5 


12.0 


1389 


Bl: V:: 


27 


70 


292.6 


06.1 


1419 


BO-0.2 III-I 


31 


98 


259.0 


10.0 


1459 


09.5 II 


34 


145 


272.7 


14.3 


1460 


B0-2 V 


34 


172 


282.7 


21.0 


1468 


09.5 V 


36 


16 


321.0 


13.3 


1500 


B0.2 III 


39 


40 


275.2 


09.0 


1531 


06 V((f)) 


43 


-25 


308.0 


09.0 


1553 


07 V 


47 


-09 


321.5 


07.3 


1563 


07.5 Il-III(f) 


47 


-04 


271.6 


06.9 


1584 


B0-1 V 


50 


-02 


320.8 


23.1 


1604 


Bl V 


55 


85 


360.2 


17.6 


1614 


05-6 V((f)) 


56 


09 


291.2 


06.6 


1618 


B0-0.2 III 


56 


128 


270.1 


07.6 


1619 


08 111(f) 


56 


102 


357.8 


20.0 


1643 


05 V 


60 


128 


279.0 


06.6 


1661 


Bl III 


62 


124 


322.7 


07.5 


1685 


BO.5-0.7 III-I 


66 


161 


291.8 


11.6 


1729 


Bl II-III 


71 


80 


283.2 


17.3 


1737 


B1.5 III 


71 


139 


339.1 


05.7 


1969 


BO. 7 IV 


113 


74 


329.6 


12.7 


1987 


B2 I 


120 


-113 


294.8 


05.0 


10001 


04 V 






246.1 


14.7 


10003 


Bl-1.5 V 






279.9 


06.9 


Suspected Binaries 


222 


09.5-BO V 


-62 


143 


198.9 


26.4 


613 


08.5 V 


-23 


-154 


203.2 


13.7 


1024 


O9-B0 V 


05 


-110 


510.7 


27.9 


1191 


BO.2-1 III-V 


13 


-30 


345.7 


28.8 


1369 


08.5 V 


26 


-09 


318.1 


78.5 


1938 


07.5 V 


107 


134 


350.0 


18.0 


1988 


B0.5 V 


121 


-22 


300.0 


18.0 



4-1- The distribution of radial velocities 

In Figure [j] we plot the spatial distribution of the mea- 
sured stars in three different velocity bins in order to check 



for the presence large scale motions, such as rotation, or 
clustering of stars in large substructures of different kine- 
matics. Unfortunately, the distribution of the points, dic- 
tated by the geometry of the spectral masks, precludes a 
finer analysis, but a visual inspection of the graph reveals 
no evidence of clumping of stars with similar velocities, 
nor of rotation along the axis defined by the observational 
technique. 

We can safely proceed, therefore, to draw histograms 
and estimate velocity dispersions. Figure |^ presents ve- 
locity histograms of the single stars in Table [j] with two 
different bin sizes, and two different origins to check for 
sampling effect. The top two panels are for a binning of 
11 km s _1 (corresponding to our estimate for the mean 
random error in the velocities) and the lower two for a 
binning of twice this error. The multiple peaks of the first 
plot disappear when the bins are shifted by 5 km s _1 (half 
a step) , indicating that they are artifacts of the small num- 
ber statistics. This is confirmed in the lower panel where 
shifting the sampling by half a step (11 km s _1 ) does not 
change the distribution in any significant way. The hy- 
pothesis of a Gaussian distribution is valid, based on % 2 
tests performed to the distribution. 

We conclude that there is no evidence for statistically 
significant peaks in the radial velocity distribution of the 
cluster. Thus, we can use all the data to estimate the ve- 
locity dispersion of the cluster. After correction for mea- 
surement (internal) errors and zero point errors (from the 
[O i] Auroral line), the radial velocity dispersion of 48 sin- 
gle stars in the cluster is 32 km s _1 . This is much larger 
that the value expected if the cluster is virialised with a 
total mass equal to the photometric mass and also much 
larger than our combined errors (<7tot ~ 15 km s _1 ). 

4-2. Binaries 

Spectroscopic binaries are very difficult to detect in a sin- 
gle observation, specially single-lined ones. Clearly, there- 
fore, the number of binary candidates listed in Table [l] 
is only a lower limit.Th.is is consistent with the fact that 
studies of young open clusters indicate 30% to be a typical 
percentage of binaries detected in these systems (Garmany 



et al. 198C, Levato et al. 1991) to be compared to 13% in 



our sample. 

Our data alone, therefore, can only be used to deter- 
mine a lower-limit to the effect of binaries in the observed 
velocity dispersion of NGC 2070. An upper limit can be 
obtained using the Montecarlo simulations of Bosch and 
Meza ( |1998| ). Assuming that all the stars in the clus- 
ter are binaries with their center of mass at rest within 
the cluster the models predict a velocity dispersion of 
Cbin ~ 35 km s _1 . This must be compared with our ob- 
served dispersion of 46.5 km s _1 including the 55 stars of 
Table [l] (corrected for observational errors). If we exclude 
star #1024 with a radial velocity of 510 km s _1 which 
may not be a member of the cluster, the dispersion is 
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Fig. 2. Distribution of radial velocities for two bin sizes of corresponding to the typical measurement error (11 km s -1 , 
top panels) and twice that error (bottom panels). The histograms are shifted by half a bin between the left and the 
right panels to illustrate features due to sampling statistics 



reduced to 36.5 km s _1 . This result is consistent, within 
the uncertainties, with the hypothesis that most of our 
observed velocity dispersion for the cluster is due to bi- 
nary motions. A (very uncertain) lower limit for the 
virial motions of the stars in the cluster potential is thus, 
cr vir = V36.5 2 - 35 2 - 10 km s" 1 . 

For comparison purposes we can estimate the 
expected velocity dispersion assuming the cluster 
is virialised and the total mass is equal to the up- 
per photometric mass limit from Paper III. From 



the density distribution derived in the same paper, 
we estimate a core radius of 0.5 pc which yields 

CTphot = 18 km s _1 . 

4-3. Mass Segregation 

In Paper II, using only the stars with spectroscopy, we 
found that the most massive stars in NGC 2070 were pref- 
erentially found closer to the center of the cluster. This was 
interpreted as tentative evidence in favor of the existence 
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of mass segregation, as was originally advocated by Malu- 
muth and Heap (1995). This conclusion was somewhat 
weakened in Paper III which presented a detailed analy- 
sis of the IMF in several rings around the cluster center. 
The IMF was found to have the Salpeter slope almost ev- 
erywhere with the exception of the very core where, com- 
bining intermediate mass HST data from Hunter et al. 
( |1995| ) with our high mass end data, we found marginal 
evidence for flattening. The most important "mass segre- 
gation" was found in a "ring" 6 pc away from the cluster 
center, again weakening the idea that closer to the center 
we would find the major relative concentration of massive 
stars. We should point out in the context that, because 
of the strong density gradient, the half-mass ratio of the 
cluster is very small. This explains the large concentra- 
tion of high mass stars in the centra l parsec of the cluster 
found by Massey and Hunter ( 1998 ). 

The two-body relaxation time for NGC 2070 [] is about 
two orders of magnitude larger than the age of the stars. 
Therefore, if mass segregation is indee d present i t mus t 
be primordial (Bonnel & Davies 1998 , Elmegreen 
In either case, dynamical or primordial, we expect to sec a 
difference in the velocity dispersion of the stars as a func- 
tion of mass in the sense of it being lower for more massive 
stars. Moreover, if mass segregation has a dynamical ori- 
gin, we expect to see energ y equ ipartition between stars 
of different masses (Spitzer [l96g| ). 

We can test these hypotheses by splitting the observed 
non-binary stars into two equal groups of 24 objects ac- 
cording to their masses as indicated by their spectral types 
(Table[|). The result is presented in Table |. 

Table 2. Mass segregation 



Mass range 


Mean mass 


Velocity dispersion 


> 23.5 M 


49.6 Mq 


27.8 km s" 1 


< 23.5 M 


19.4 Mq 


36.7 km s" 1 



The Fischer F-test on both distributions gives a value 
of F=1.6 corresponding to a probability of 27% that 
both samples are drawn from the same parent distribu- 
tion. Thus, there is tentative, but not conclusive, evidence 
that the massive stars have a lower dispersion. The ra- 
tio of mean energy (M 2 ) between the two mass bins is 
~ 1.5 ± 0.1, significantly different from the equipartition 
ratio, r = 1. So if the mass segregation is indeed present, 
it is most likely not due to two-body relaxation. We re- 
mark, however, that our radial velocity data samples very 
sparsely the inner 10 pc of the cluster, where we concen- 
trated our photometric study, and which contains most of 
the cluster mass. 



5. Conclusions 

In spite of the relative low spectral resolution, our data al- 
ready provide important new results about the dynamical 
state of the 30 Dor cluster. First, the velocity dispersion 
is much too large to be due random motions of the stars 
in the gravitational potential of the cluster. Instead, the 
observed dispersion seems to be entirely dominated by bi- 
nary orbital motions. Thus, the first important results is 
that higher spectral resolution alone is not sufficient to 
probe the dynamics of the cluster; it is also crucial to 
have observations with good time resolution in order to 
find (and exclude) binaries. Second, there is no strong ev- 
idence for dynamical mass segregation in the sense of mas- 
sive stars moving with lower random velocities. If present, 
the effect is masked by binaries, so again, it is crucial to 
obtain data for several epochs. Finally, the virial dynam- 
ical mass of the cluster is comparable within factors of a 
few with the photometric mass. Therefore, using a reason- 
able number of non-binary stars it should be possible to 
place useful constrains on the IMF slope below 1 M Q . 

The strong conclusion of this investigation, therefore, 
is that it would be very worthwhile to obtain time resolved 
spectroscopy of a sample of 100-200 stars in the cluster. 
The FLAMES integral field spectrograph on the VLT ap- 
pears ideally suited for such study. 
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